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Abstract. In this paper it is shown that the steady-state weights of the asymmetric 
simple exclusion process (ASEP) with open boundaries and parallel update can be 
written as a product of a scalar pair-factorized and a matrix-product state. This type 
of state is also obtained for an ASEP on a ring in which particles can move one or 
two sites. The dynamics leads to the formation of an excess hole that plays the role 
of a defect. We expect the process to play a similar role for parallel dynamics as the 
well-known ASEP with a single defect-particle (that is obtained in the continuous- 
time limit) especially for the study of shocks. The process exhibits a first-order phase 
transition between two phases with different defect velocities. These are calculated 
exactly from the process-generating function. 
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1. Introduction 

The asymmetric simple exclusion process (ASEP) has been used to model different 
dynamical systems such as traffic flow and biological processes. It is originally defined 
in continuous time on a discrete one- dimensional lattice. Particles on the lattice can 
move one site to the right at a certain rate if the target site is empty (see e.g. [lj for 
a review). The model with periodic boundary conditions is known to have a uniform 
stationary measure pD|. However employing open boundary conditions (particles enter 
the system at the one end and leave it at the other end of the lattice at certain rates) 
leads to so-called boundary-induced phase transitions. The steady state is of the matrix- 
product form [Bj. There are, however, some generalizations of the ASEP with periodic 
boundary conditions with phase transitions. An example is the ASEP with a single 
defect particle that can itself move forward on empty sites and can be overtaken by 
normal particles [T2| H]. 

The ASEP has been extended to various discrete updating schemes rather than a 
random-sequential update [U [3], the probably most important of which is the parallel 
update which is typically used for traffic flow simulations [2] . Parallel means that all the 
sites are updated simultaneously and particles attempt a hop forward with probability p. 
The introduction of such a parameter is necessary to interpolate stochastically between 
purely deterministic movement (p = 1) and the continuous-time limit (p — > 0). For open 
boundary conditions the parallel-update ASEP could be solved by two sophisticated 
versions of the matrix-product ansatz [SI IE], see section 2. However since then there is 
somehow a lack of new exact solutions for steady states of cellular automata. 

The outline of the paper is as follows. Section 2 shortly reviews the exact solution 
for open boundaries as obtained in [S].Then we give an alternative form of the steady 
state that is a product of a scalar pair-factorized and a matrix-product state which is 
the most straightforward form one would expect from the knowledge of the solution for 
periodic boundaries and the open-boundary solutions in other update versions. It is 
shown how it corresponds to the solution in [8]. In section 3 we investigate a process 
on a ring in which particles have maximum velocity two and calculate its exact steady 
state. In section 4 we perform an alternative type of proof as in the previous papers 
[SI E] ■ Then, in section 5 we calculate the normalization-generating function and extract 
the phase behavior and phase-typical asymptotic quantities are compared to computer 
simulations. Finally in section 6 we map the process in the continuous-time case onto 
the defect ASEP [121 [TJ]. 

2. The asymmetric simple exclusion process with parallel dynamics 

The ASEP is defined on a one-dimensional lattice with L sites, enumerated I = 
1,2, ... ,L. Each site I may be in one of two possible states (expressed by a local 
state variable r), namely either occupied by one particle {ji = 1) or empty {ji = 0). 
In the case of open boundaries the system is coupled to two boundary reservoirs, 
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one to the left and one to the right. A particle enters onto the first site if it is empty 
with probability a and a particle may exit from site L with probability j3. In the bulk 
a particle can move one site to the right with probability p if the target site is empty. 
Note that every site is updated simultaneously. The exact form of the stationary state 
can be obtained by the matrix-product ansatz originally introduced for continuous time 

151 ■ 



2.1. The site-oriented ansatz 

By introducing boundary vectors (W\, \ V) and matrices E and D for holes and particles 



respectively Evans, Rajewsky and Speer [5] could show that 

P(n, ...,r L ) = Z L \W\ n [(1 - n)E + r t D] \V) (1) 

i=i 

gives the correct steady state when the operators satisfy the bulk relations 

DDEE = (1 - p)DDE + (1 - p)DEE + p(l - p)DE, (2) 

DDED = DDD + (1 — p)DED + pDD, (3) 

EDEE = (1 — p)EDE + EEE + pEE, (4) 

EDED = EDD + EED + pED, (5) 

as well as relations for the right boundary 

DDE\V) = (1 — p)DE\V) + DD\V), (6) 

EDE\V) = ED\V) + EE\V), (7) 

DD\V) =^^D\V), (8) 

ED\V) =ZE\V), (9) 
and left boundary 

(W\DEE = (1 — p)(W\DE + (W\EE, (10) 

(W\DED = (W\DD + (1 — p)(W\ED, (11) 

(W\EE = P ^~ a \ w\E, (12) 

a 

(W\ED = ^(W\D. (13) 



Note that - as a consequence of the particle-hole symmetry of the process - this relations 
are invariant under exchanging a <-» /3, E «-> D, (W\ <-> \V) and at the same time 
inverting the order of the enumeration of cells (site % is replaced by site L — i + 1). 
The ansatz (CQ) together with (121 T131) are a notation for certain recursion relations for the 
weights on system size and particle number [5]. The authors wrote for simplification of 
the calculations the ansatz 



E = 




gDi \ 










o ) 







(14) 
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(w\ = «wi|, 7 - ?L T - mi \V) = \\Vx),-^-\V x )), 

with a certain constant g > 0. Here the matrices E and D are effectively rank four 
tensors since its components are itself matrices. It turns out that the given operators 
flUD fulfill the quartic algebra O H3]) if 

D 1 E 1 = (l-p)[D 1 + E 1 +p], (15) 

(Wi\E, = ^"^ (^il, (16) 

D^^-^W. (17) 

All physical quantities can be expressed through D 1} Ei and (Wi\, \Vi). However using 
(fj~4l) the weights (DD) become difficult expressions, namely a complex sum over matrix 
elements, that do not have an obvious physical meaning. Therefore it would be desirable 
to have a more easy formulation of the weights. 

2.2. Alternative solution (Product of a scalar pair- factorized and a site-oriented 
matrix-product state) 

Before we present an alternative formulation of the weights, let us remember the type 
of solution for periodic boundary conditions. The time evolution on a ring leads to a 
pair-factorized stationary state with weights [7]: 

Fr(r 1 ,...,r L ) = f[t(r l ,r l+1 ) (18) 
i=i 

with some simple two-site factors t(rj,Tj + i) obeying 

£(01) V yj t(00) V ; 

A useful choice is 

£(11) = (1 - p)t(01) and t(10) = t(00). (20) 
We rewrite f fTSl matrix product state 

Fr s (r 1; . . . , ti) = tr n [riD + (1 - n)E] (21) 
i=i 

Of course 2x2 matrices D and E of the form (fill) solve also the periodic case [5J [18] but 
we prefer to write the matrices in an ordinary vector basis, the t(ra) becoming matrix 
elements (t\(E + D)\a) in the style of an Ising transfer matrix: 

V t(io) o j ' v ) 

which together with (12"T|) is the simple matrix equivalence to ([TBI . In this representation 
it is obvious that the product is self-consistent, i.e. that no terms . . .t{Ta)t{r'a') with 
a ^ t' occur and leads to one single term (IT8|) . 
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Now we come back to the case of open boundaries. Inspired by (1221) we take 
alternatively to ([Tj 



_ / t(00)Ei \ / f (01) A \ 

^-U(io)^ oj' D -(o t(ii)A J- (23) 

Here t{ra) are the two-site factors of the solution for periodic conditions ( 120|) . We set 
£(00) = t (10) = 1 according to (1201) . Since the operators E or D of the form (1231) have 
the structure of (1221) with factors £(tj_iTz) in a matrix for site /, the correct connection 
to the operator for site I — 1 is guaranteed. The boundary vectors (W\ and \V) read 

(W\ = Wil—mi, (24) 
p — a 1 — p 

and for the bulk we find t{Ql)D x Ei = t{ll)D l + E x + p. So setting £(01) = (1 - p)- 1 
and £(11) = 1, recovers (fT5lh Using the new operators fl22|) it is rather obvious that the 
ansatz (CQ) yields formally 

L 

F(n, . . . ,r L ) = w(ri)t(ri,r 2 ) . . . t{r L -x, r L )v(r L ) x JJ [ n D x + (1 - r 2 )^i])K) (25) 

i.e. a superposition of a pair-factorized state (reflecting the nearest-neighbor correlations 
of the parallel update) and a matrix state (as for other discrete-time updates such as 
ordered sequential and sublattice-parallel updates [HE]). Here t(ra) is defined through 
(|20|) and the boundary factors are 

w(n) = t(0ir + -^t(ll) T \ v(tl) = \ Z ^ L - (26) 
p — a I - (Jt l 

The pair-factorized pre-factor obviously distinguishes the states between boundary-site 

occupations and (as on the ring) by the number of 01-pairs (or 11-pairs respectively). 

2.3. Connection with the bond-oriented ansatz 

De Gier and Nienhuis [8] alternatively solved the parallel ASEP with open boundaries 
through a bond-oriented matrix ansatz: 

P(n, r 2 , ... , t l ) = (W(ri)|M(n, r 2 )M(r 2 , r 3 ) . . . M(t l ^t l )\V(t l )) (27) 

The vectors and matrices M(to~) are 

M( s ( (1 - r)(l - a)M{00) (1 - r)^(01) \ 

M (^) = ( r (l- (7)^(10) r^(ll) J (28) 

(H/(r)| = ((l-r)(W(0)|,r(W(l)|), (29) 

|V(r)) = ((l-r)|V(0)),r|V(l))) i . (30) 

We now give a relation between the site-oriented and bond-oriented solutions. In (J23J) 
we have to take £oi = 1 and tii = 1 — p. Then the connection is: 

A4(00) = M(10) = E h M(ll) = (1 - p)M(01) = Di (31) 

(W{0)\ = (W 1 \E 1 + (W 2 \E 1 , (W(l)\ = (W 1 \D 1 + (l-p)(W 2 \D 1 , (32) 
|V(0)> = |Vi>, |V(1)> = \V 2 ), (33) 
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and therefore 

E = M(00) + M(10), D = M(10) + M(ll), (34) 

((w|£?,o) = (w(o)|, (o,{^p) = TOI> |v> = |v(o)> + |y(i)>. (35) 

In fact one sees that our choice ( |22l) is closely related to the bond-oriented solution and 
is just rewritten systematically in terms of a site-oriented matrix product. 

3. Solution of an ASEP on a ring with excess-mass formation 

We are going to consider an asymmetric exclusion process on a ring with L sites, N 
particles and periodic boundary conditions (site L + 1 = site 1). The system evolves 
under parallel dynamics according to the local update rules 

100 -> 001, with probability p, (36) 

101 — > 011, with probability (3. 

It turns out that the stationary distribution of probabilities for the possible config- 
urations is not ergodic, i.e. only a subspace of configurations is reached as the time 
increases. This stationary distribution depends strongly on the parity of the number 
L — N of holes (unoccupied sites). We specify a certain configuration of particles by 
the set of gaps (number of holes) between the particles: C = {rii, n 2 , ■ ■ ■ , n N ). The 
model dynamics is such that odd-valued gaps can not be constructed, however they 
can turn into even gaps when a configuration C(. . . l[any odd number of 0s] 101 1 . . .) 
moves with conditional probability (3 into a configuration with two odd-valued gaps 
less. These processes appear again and again until there remain either no more odd 
gaps (L — N even) or exactly one odd gap (L — N odd). In the latter case there remains 
so to speak a single excess hole (comparable to the concept of excess mass in the math- 
ematical literature). In the following we are going to consider these two cases separately. 

For even number of holes the system arranges such that only even-length gaps 
remain. The weight for a configuration factorizes into N factors, one for each gap. All 
positive even gaps have the same weight. Only the weight for zero gap is different: 

N 

F(ni,n 2 ,...,n N ) = Y[f(n tt ), (37) 

with 

{1 — p, for n — 0, 
1, torn = 2,4,..., (38) 

0, for n — 1, 3, . . . 

In the subspace of even gaps this is equivalent to the solution of the usual ASEP on a 
ring (TIB"]) which can simply be written as (|37j) . with /(0) = £(11) and f{n) = £(01), for 
n > 1 being a possible choice j9] . So this case is simple and well-known and will not be 
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discussed here further. 



For odd number of holes there remains only one odd- valued gap (Configurations with 
more than one odd- valued gap have probability zero in the steady state). We introduce 
a matrix-product ansatz for the weight of particle p, being followed by 2n fJL + 5^n holes, 
H — 1, . . . , N. In contrast to the usual formulation where a matrix E represents a hole, 
the matrix E here stands for a pair of neighboring holes (00). D represents particles (1) 
and |V)(W| = A stands for the excess hole together with the particle to its right (01). 
So we use here the convention that the position of the excess hole is always at the right 
end of the gap (00 00 ... 00 1) between particles N — 1 and N. The ansatz reads: 



F(2n 1 ,2n 2 ,...,2n N + l) = (W\ 



N-l 

Y[ E n »D 
n=i 



E nN \V). (39) 



We note that it is possible [10J to write a matrix ansatz allowing for any number of odd 
gaps that leads finally to ( J39l) . However from our argumentation above it should be 
clear that only one odd gap remains and so we base our solution on this simple finding. 
From diagonalizing small systems we find a quartic algebra for the process related to 
(l2lfT3l). In fact by transforming the matrices E — > (1 — p)~ 1 E and D — > (3D almost the 
complete set of relations can be mapped onto the algebra of the parallel ASEP with open 
boundaries. To be precise one recovers exactly (|2lfT2l) with a — p, the only exception 
being ([TBI which has to be replaced by 

(W\ED = (W\(D + p). (40) 

This is in accordance with the dynamical rules ( 136|) leading to the fact that even for 
p = (3 the particle-hole symmetry is broken. The matrix transformation for D and 
E mentioned above can be omitted in the calculation since for fixed particle and site 
number it leads only to an overall factor in the normalization constant 



E ■■■ E %>„,mW 



J N,M - ■ ■ ■ °][> M ,M ( 

ni=0 njv=0 



N-l 



n En " D 



E nN \V) (41) 



for TV particles and 2M + 1 holes according to ( |39l) . This factor is well-defined if p ^ 1 
and (3 ^ 0. So we shall omit these in fact less interesting cases here. The initial values 
are 

(W\E\V) = ( 7 + P)(W\V), (W\D\V) = py/P(W\V), (42) 

for some constant 7 > 0. We just note that this leads to (W\DE\V) — (1 — 
(3){W\D\V) + {l-p){W\E\V)+p'Y{W\V}. We found that the algebra (JKUIIiD can not 
be simplified by ( |T4"1) . However with the alternative ansatz ( 1231) it can. So we expect 
that (1231 is the more robust form that holds even for broken particle-hole symmetry in 
probabilistic cellular automata that give rise to a matrix-product state. We take again 
t(00) = £(10) = t(ll) = 1 and £(01) = (1 — p)~ x . Then the weights can again be written 
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as a superposition of a pair-factorized and a matrix state as (1251) . We can rewrite 



as 



F{2n u 2n 2 , . . . ,2n N + 1) = (W\ 



N-l 
Ml =1 



E nN \V) 



1 -p5(njv,0) 



.(43) 



I - P5(n N , 0) 

In our opinion this form of the weights helps to understand the solution of this type of 



models. However it is not less convenient to work directly with the matrices (1231) which 
read here 

' {l-p)- 1 D l \ 
,0 D x 



E 



and lead to boundary factors 



D 



(44) 



(W\ = (0,(W 1 \), \V)=[\V 1 ) : 



1-P 

1-/3 



This choice leads here to a ternary algebra for the indexed matrices: 

D 1 E 1 =(l-p)[D 1 + E 1 +p], 
p(l-P) 



E 1 D 1 \V 1 ) 

D 1 D 1 \V 1 ) 
(\Vt\EiEi 



p{l-/3) 







DM), 
l-p)(W 1 \E 1 , 

l-p)(W 1 \D 1 +p(l-p)(W 1 



(W 1 \E 1 D 1 

Translating (j42p into the form with indexed matrices gives (W / i|£'i|V / i) 
p)/{\ - 0) {WM) and (Wi\D x \V x ) = pn/P (WilFi). A useful choice for 7 (which 
coincides with the representation fl56|) that we give below) is 7 = 1 — (3. Then one has: 



(45) 

(46) 
(47) 

(48) 

(49) 
(50) 

(7 + /*)(! - 



<Wi|£i|K> 



1 — p 



(WxlVx) and (Wi|Di|Vi) 



p{\-(3) 







The choice 7 = 1 — (3 is useful because the algebra (J4b|) simplifies to 

=(l-p)[D 1 + S 1 +p], 



Al^i) 



p(l-ff) 



<Wi|Vi). (51) 

(52) 

(53) 

(54) 
(55) 



{l-p)(W x \E u 
(W 1 \E 1 D 1 = (1 -p)(W 1 \D 1 +p(l -p)(W 1 \ 

since here the first two rules are quadratic (and are the same as for the open-boundary 
case) . 

For these relations we find the representation: 

/ 

(1-p) 


0... 



£1 



V 



(56) 



/ 
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D 1 



( p(l 




V ••• 



P)/P p/(3 

(1-p) 







p/p 

1 

(1-p) 





p/p 

1 
1 

(1-p) 



(W x \ = (1-/3,1,1,1,...), |^) = (1,0,0,...)' 



(57) 



(58) 



It should be mentioned that the known representations for the open-boundary case [8] 
(which reflect the particle-hole symmetry by a symmetry in Di and E\ as well as (Wi\ 
and |Vi)) can not be generalized to represent the present process. However f )56ll58l) can 
be changed to represent the open-boundary case, namely by changing (Wi|l) from 1 — f3 
to 1. 



4. Proof of the steady state 

Different techniques have been used previously to prove the stationary states of the 
ASEP with parallel update and open boundaries. The so-called canceling mechanism 
could be generalized for parallel dynamics [I], however it may involve more than two 
neighboring sites [8]. Here the problem remains to find representations for auxiliary 
matrices as well which is a difficult task in general. In [5] the state was proven by 
using the quartic algebra in a rather more mathematical language. We will prove the 
ansatz here in an alternative way, namely by using the quadratic and cubic rules for the 
matrices D\, E\ and A\ instead of the quartic rules for D, E and A. To do this we will 
derive the master equation from the local dynamical rules and afterwards will prove it 
by using the rules for DiEi, DiAi, A\E\Ei and A\E]_D\. 



4-1. Derivation of the master equation 



It is not obvious how to write down the master equation here. We now write the state 
of the system as the ket- vector |n 1; n 2 , ■ ■ ■ , n N ), denoting particle 1 followed by n\ holes 
and so on. This may formally be obtained by the tensor product of the single-particle 
states \rifj). Let djk(n^) be the transition probability for particle // to go from state 
\ n fj, + j + k) into \nn) on moving j cells while particle \l + 1 moves k cells. Then the 
master equation can be written as (compare [9]) 



N 



(F\{n,}) = (F\trl[Tl 



n 



with the transfer matrix 

/ doo(^)l^) doi(n M )|7i M -l) 

^io(^ M )l^ + 1 ) ^n(^)l^) 
d 2 oK0l?V + 2) d 2 i(n lx )\n lx + l) 

The transition probabilities follow from (1361) and are 

d ok (n) = 4,n + (1 - P)6 k , n -i + (1 - p)9(n - 1 - k) 



d 2(n^)\n fl - 2) ^ 
di2(n^)\n^ - 1) 



(59) 



(60) 



(61) 
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dik(n) = P6 k>n , (62) 
d 2k (n) =p6(n- k + 1). (63) 

Since we know that in the steady state there remains only one odd gap between 
the particles we use this to simplify the equation. Let the odd gap be between 
particle N and particle 1. Then we ask for the probability flow into the state 
|2ni, 2n 2 , . . . , 2njv-i, 2n N + 1). To obtain this state either particle N or particle 1 have 
been in the odd state before, since the odd gap can move only backwards. All other 
particles have been in an even state. Using this one finds for T(2n-\): 

/ (S Jllfi +(l-p)9(n 1 ))\2n 1 ) (<S ni)1 + (1 - p)9(n 1 - l))|2m - 2) \ 
T(2m)= pS nu0 \l) P6 nul \l) 

v p|2ni + 2) p0(m)|2ni) 

The second column vanishes because particle 2 can not have moved one site since it had 
an even gap in front as claimed before. Now using the matrix ansatz this can be written 

as 



(64) 



/ 



T(2m) = 

Equivalently one has for the bulk 



ni+l Di 









n 1 )E^~ 1 D, 

r 



(65) 



T(2n fl ) 



( E^D X 




V 



— r> 1 1 







H = 2,...,N-1, 



(66) 



and for T(2un + 1): 

/ (l-p^Mj 



T(2n w + 1) 









^_Er +1 D 1 P 6{n N )ErA l ) 



(67) 



Note that the component of the second row and second column, containing a factor d n , 
vanishes in every transfer matrix since it is impossible that a particle and the particle 
in front of it move at the same time only a single site in the steady state. 

Now inserting these matrices into the master equation one ends up with a product 
of bulk transfer matrices that is rather difficult to handle. The crucial step in deriving 
the master equation is the following similarity transform: Take 



/ (l-p)E l 
L = 

V -pEi 

Then one has LR = RL — E 1 G 
( E n 1 » +1 D 1 + ^ 









(1 



(1 



p) \ 

p) ) 



and R 



( 1 










-1 \ 







(68) 



LT(2n^)R 




V0 



5 1 and the convenient expression (for \x = 2 . . . N — 1): 

ein^E^D.E, -E^ +1 D 1 + 6(n lx )E^D 1 E 1 \ 


J 



.(69) 
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From here a straightforward calculation involving successive simplifications (without 
using the algebra (I52H55I) ) shows that the master equation can finally be written as 



7V-1 



"N 



N-l 

[1 - /3<W, " P°M] trA: J] [(1 " p9(n fl ))E^D x + p9(n^E^- 1 D x E x ] E\ 



pnjv 



1 -p9(n N ) 
1 - /30(n*) 



JV-1 

trAi (5 ni)0 + p5 ni)X E x ) II [{l-pe{n tM ))E^D l +p9{n^)E^- 1 D 1 E 1 \ E X N D x 

fj,=2 



N-l 



+ P 6{n N )trA 1 E 1 J[ [(1 - P e{n ll ))E%"D 1 + p9(n, jb )E^ ^ 



(70) 



Note that only the simple structure of (1231) allowed for a closed expression of the master 
equation in terms of the primed operators. 



4-2. Proof of the matrix-product ansatz 

In the following we assume always N >2, since the case N = 1 is trivial. For the proof 
the following simplification of the bulk terms under the product is essential: 

(1 - p9(n fi ))E^D 1 + p9(n„)E^' l D x Ei = 5 nfl>0 D x + 9(n„)E^~ l [(1 - p)E 1 D 1 + pD x E x ] 
= 5 nfl , D x + (1 - p)9{n ll )E^- 1 [E X D X + p(D 1 + E x + p)\ 



5 n ,D x + (1 - p)9(n„)E^ 1 [(Ei + p)D x + p(E x + p)\ 



= 5 n ^D x + (l-p)9(n tl )E^- 1 (E l +p)(D l +p). 
Here we have factors (E x +p)(D x +p). Note that from ( 1521) it follows [5] that 
(1 + +p) = D X E X . 



(71) 



(72) 



This can be used to simplify the following equation which turns out to be the key to 
the proof: 



(D x +p) 5 n ^D x + (l-p)9(n fl )E^- 1 (E x +p)(D x +p) 



5 n nD x + (1 -p)9(n„)(D x +p)E^-\E 1 + p) (D x +p) 



= KfiDi + (1 -p)9(n lM )(D x +p)(E 1 +p)E%" '\ (D x +p) 

= 'KfiDi + 9(n^D x E x »] (D x +p) 
= D x E^(D x +p). 

Here we have used the fairly simple but essential commutation relations E x (E x + p) — 
(E x +p)Ei, Di(Di + p) = (Di + p)D x . As a consequence one has 

(D x +p)H[(l-p9(n l ,))E^D x +p9(n l ,)E^- 1 D x E x \ = JJ [D X E^] (D x +p). (74) 
fj, /i 

In the following we consider only the case un = since the case > can be handled 
in a similar fashion [10]. For = we simplify the master equation ( J70l) on both sides 



(73) 
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N-l N-l 

ptxAiE? J] [^1^1 =/3trAx J] [(l- P 9(n^)E^D 1 +p9(n^)E^- 1 D 1 E 1 
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JV-l 



+p(3trA 1 (8 nu0 +p6 nul E 1 ) ]J [(1 - p9(n fl ))E^D 1 +p9 1 D 1 E 1 



H=2 



(75) 



In the first term on the right-hand side (rhs) the factor corresponding to \i = 1 is 
extracted from the product. Rewriting it with the help of (ITT]) and combining terms on 
the rhs with rti = yields 



7V-1 



. . . = p5 nifi trA 1 (D 1 +p) J] [(1 -p9{n fl ))E^D l + p9(n^)E^ X D X E X 

N—l 

+P(l-p)9(n 1 )trA 1 E^-\E 1 +p)(D 1 +p) [(1 - P 9(n^)E^ D x + p9{n ll )E^- 1 D 1 E 1 



N-l 



+ P 2 p6 m>1 txA 1 E 1 J] [(l-p9(nJ)E^D 1 +p9(n^)E^ X D X E A 



y,=2 



(76) 



In the same way one can extract the factor for \i = 2 in the third term of (1761) and rewrit- 
ing it with dTTJ). The resulting factor A X E X [S n2 ,oDx + (1 - p)9(n 2 )E^ 1 {E 1 + p){D 1 + p) 
is for n 2 = 0, due to (|53|) . equal to (1 — p)A 1 (D 1 + p). For n 2 > it becomes (1 — 
p)AiEi 2 (Di+p) as a consequence of (15^1) . Concluding one finds after combining both ex- 
pressions A X E X [8 n2fi D x + (1 - p)9(n 2 )E^-\E l + p)(D 1 + p)] = (l-p)A l E^(D 1 +p). 
Then (1761) gives 

N-l 

... = i38 nifi tiA 1 {D 1 +p) J] [{l-p9{n fl ))E^D 1 + P 9{n fl )E^- 1 D 1 E 1 



H=2 



N-l 



+ 



(3(l-p)9(n 1 )tiA 1 E^- 1 (E 1 +p)(D 1 +p) [(1 - P 9{n^))E^D 1 + P e(n fl )^ 1 D 1 E 1 



fi=2 



N-l 



+ P 2 (1 



/i=3 



(77) 



Now use f!7^|) and /3(.Di + p)Ai = p^i which follows from (153]) and combine terms with 
ni = 1. Then the master equation turns into 

N-l N-l 

ptiA.E^ J] [DiE^] = pS^otrAx ]J [d^ 



y,=2 



li=2 



N-l 



+p(l -p)S nul tiA 1 [(E 1 +p)D 1 +p 2 ] E? 1] [ D ^ E i" 

/i=3 
N-l 

+p\l-p)9{n 1 -l)tiA 1 E^-\E 1 +p) J] [D x El 



fj,=2 



(7f 



In the second term on the rhs there is a factor (1 — p)Ai [(Ei + p)Di + p 2 ] which can be 
simplified with (155]) and yields simply A\E\Di as one can check in a few lines. Finally 
consider in the third term on the rhs (1 — p)9(ni — VjAiE'l 1 " 1 (E 1 +p). Since this term 
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only for n\ > 2 gives non-vanishing contributions, (I54p can be applied and leads to 
(1 — p)A 1 E™' 1 l {Ei + p) = AiE^ 1 . Inserting these results in ( 1751) yields the required 
identity. 



5. Asymptotic behavior and phase transition 

As claimed above the process with even number of holes corresponds to the usual ASEP 
which is well studied so again we focus only on the case with a single excess hole. In 
contrast to the open boundary ASEP, on a ring we have a fixed number of particles and 
holes. The calculation is done grand-canonically by introducing fugacities x and y for 
particles and hole-pairs respectively. Consider the grand-canonical probability p_ (n) of 
finding a particle directly behind the 01-pair while there are a total number of n other 
particles and hole pairs: 

M»> = I<tl/|C "" g|t/> (79) 

where 

C = C(x,,)=xD+yE=(** ^"f'M ( 80 ) 

\ yEi xD x J 

results from (1231) . We note that this can be related to the corresponding expression that 
one would obtain from (jUj) by a simple similarity transform [10]. The grand-canonical 
normalization for an excess-hole system of a total number of n + 1 particles and hole 
pairs is Z n = (W\C n \V). The nominator in (179]) can be simplified: 

(W\C n D\V) = (WlC^xDDlV) + (WlC^yEDlV) 

- p(1 ~ ^(WlC^xDlV) + ?-{W\C n - l yE\V) 



= j-(W\C n \V) - px{W\C n - 1 D\V). (81) 

We define S n = (3/p (W\C n D\V) , so that Z n = S n +pxS n -i, for n > 1. The asymptotic 
form of S n is always S n ~ A - ™ (which follows from the theory of generating functions, 
see [II]), with a site-representing fugacity A, so that Z n ~ (1 +px\)\~ n and 

pxX 

p - = p(i+ P x\y (82) 

Explicit derivation of the generating function S = J2 S n X n and analyzing its singularities 
shows the existence of two phases (see Appendix A). 

• First phase 

The first singularity results from a pole. One finds a relation between the fugacities 
in the form x(X): 

8 3-p + p(l-p)X 
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Second phase 

There is also a square-root singularity in the expression for S leading to 

2 



W (1 -p)/A( 1+ pA) . 
{ 1-P(1+PA) ) ■ m 

We introduce a formal asymptotic density p ~ N/(N + M) for a system with 
iV particles and M hole pairs in the normal ASEP picture, where each matrix E 
represents a single hole. Results can easier be expressed in this form and symmetries 
as well as comparison with known results are more obvious. Instead of having a 
relation A(x), we have x(X) which fixes the density p: 

M - (85) 



Using this and equating relations (j83|) and fl84l) leads to an expression for the critical 
density: 

1(1 - ) 
For p_ we find in phase 1: 

(l-p)(l- Vl-4pp(l-p)) 
P ~ (P - p)(l - ^l-4pp(l-p)) + 2p(l - /3)(1 - p) ' 
In phase 2 the result in terms of p is tedious. Parameterized in A it reads 

P / l-V(l-p)/A(l+ P A) V 
P " /3(1+PA) ^ ^p)JX-p ) ' 1 J 

As an example we take p and (3 such that p c = 1/2, namely p = 3/4 and j3 = 1/2. 
In figure [1] one sees how the curves corresponding to the two phases fit together to the 
dotted curve coming from a computer simulation with L = 1000. For p < 1/2 the 
system is in phase 2 and for p > 1/2 it is in phase 1. One sees that the exact solution 
for L — > oo has a kink at p = 1/2. So its derivative there has a discontinuity. 

For p = /3 = 3/4 the system is completely in phase 2. The comparison between computer 
simulation and exact solution is shown in figure [3 

The question may arise if there is a choice of p and /3 for which the system is in phase 
1 for all p. However considering (|86|) shows that for the allowed parameter values p ^ 1 
and /3 ^ this is not possible. 

Now we come to the probability for the occupation in front of the 01-pair: 

_ x(W\DC n - l \V) _ ! _ yjWlEC^lV) 
P+[ } (W\C n \V) (W\C n \V) ' 1 } 
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Figure 1. The solid curves show the two solutions for p-{p) in the two phases in 
comparison with the dotted curve coming from a computer simulation. The hopping 
probabilities are p — 3/4 and [3 = 1/2, so that the critical density is p c = 1/2. The 
system size is L = 1000. 




Figure 2. p-(p) for p = (3 — 3/4 and system size L = 1000. Dotted curve: computer 
simulation, solid curve: solution of phase 2. 



Start with the nominator: 

(W\EC n \V) = (W\ExDC n - l \V) + (W\EyEC n ~ l \V) 

= x(W\{D + p)C n ~ l \V) + y{l -p)(W\EC n - 1 \V) 
= (W\C n \V) +px(W\C n - 1 \V) -py{W\EC n - l \V). 
Now define T n := (W\EC n \V) . Then one has 
T n + pyT n - X = Z n + pxZ n -i 

and we conclude that for n large T n scales as A~ n • (1 + px\) 2 / (1 + pX). Thus 

1 + px\ 



p+ = l-y\ 



1 + pyX 



(90) 
(91) 

(92) 
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This leads in phase 1 to 



p-P 2p(l-p)-l + Jl-ipp(l-p) 
P+ = ~YT\ a\ / = ( 93 ) 



and in phase 2 in terms of A simply: 

1 -p{l+p\) 



= . (94) 



O.S 




Figure 3. The solid curves show the two solutions for p+{p) in the two phases in 
comparison with the dotted curve coming from a computer simulation for p = 3/4 and 
f3 = 1/2, so that the critical density is p c = 1/2. The system size is L = 1000. 

Comparing the relations for p_ and 1 — p + one sees that due to the broken particle-hole 
symmetry there is no proper symmetry between the two relations. 
The velocity of the defect in the two phases can be obtained by 

u=p(l-p + )(l-/3p_)-/3p_. (95) 

To be precise, this is the velocity of the single excess hole. If it has a particle directly 
behind it jumps backwards with probability beta which leads to the second contribution 
—f3p-. If it has no particle in front (probability (1 — p + )) it can jump forward with 
probability p unless it also has a particle behind which moves forward with probability 
(3. This leads to the first contribution p(l — p+)(l — /3p-). Note that we always argue 
in terms of a density p that treats the hole pairs as single holes. 
Using fl82|) and fl92l) gives rise to the following expression for v: 

= py^ p^x = pKv - x ) / 96 n 

1 + pyX 1 + pxX (1 + pxX)(l + pyX) 
One sees that, due to the symmetry in x and y, the average defect velocity vanishes for 
equal densities of particles and hole pairs (v(p = 1/2) = 0). In phase 1 fl96|) is rewritten 

as 

p(p-f3)(l-p)- p (l-{3)J 



v(p) 



p(l-P)(l- p )-(p-f3)J 



(97) 
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where J is the total particle current 



J(p) 



1 



4pp(l - p) 



(9* 



which is expected since the flow should equal the result for even total number of holes 
and it has to be phase independent in our process. Note that in phase 2 the results for 
P-, p + and v are independent of (3. Figure 0] shows the exact defect velocity v(p) for 
p = 3/4 and two different values of (3. For p < 1/2 the velocity is independent of (3 
and the system is in phase 2. For p > 1/2 and (3 = 3/4 (= p) the system remains in 
phase 2 (lower curve). For p > 1/2 and (3 = 1/2 the system is in phase 1 (upper curve). 
At the critical density p = 1/2 there is a discontinuity in dv/dp indicating a first-order 
transition. This is expected, since the model for random-sequential dynamics shows the 
same type of transition [17]. This model can itself be mapped onto the ASEP with a 
defect as is explained in the next section. From this mapping also the physical reasoning 
of the phase transition should become clear. 




Figure 4. Exact defect velocity for p = 3/4 and two values of (3. The regime p < 1/2 
is independent of (3. For p > 1/2: upper curve (3 = 1/2 and lower curve (3 = 3/4. 



6. Continuous-time limit, connection with the defect ASEP and shock 



profiles 



When the hopping probabilities are so small that on average only one particle moves 
during a time-step the parallel update turns into the random-sequential update which 
mimics continuous time. To be precise, one has to replace (3 — > (3dt, p — > pdt and 
afterwards taking the limit dt —>■ 0. Then (|36l) turns into 100 — * 001 at rate 1 and 
101 — > 011 at rate (3. The algebra ( 15211551) becomes equivalent to the DEHP algebra [6] 
(here with a = 1): de = d+e, (w\e = (w\, d\v) = /3 -1 |t>). The non-vanishing weights 
can simply be written as [T7j F(2ni, . . . , 2tin + 1) = (w\ n^=L 1 ( e?Vc O eriJV 1^)- This in turn 
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is the steady state of the ASEP with a single defect particle [12] for a = 1. So how 
comes this along? The defect ASEP is defined by the local transitions: 10 — > 01 at rate 
1 and for the defect particle 2: 20 — > 02 at rate a = 1. Normal particles can overtake 
the defect: 12 — > 21 at rate (3. In fact our process can for continuous time be mapped 
onto the defect ASEP. An arbitrary stationary configuration has exactly one excess hole. 
Remember our convention that this excess hole is always localized at the right end of 
the cluster of holes to which it belongs. The mapping is as the matrix ansatz suggests: 
the 01 pair is the defect 2, the other hole pairs 00 become single holes and the other 
particles remain normal particles [T7]. Note that the particle to the right of the excess 
hole changes with time. However since the particles are indistinguishable this has no 
effect. 

The defect ASEP was first introduced and studied for a — (3 = 1 in which it 
is referred to as a second-class particle (see [13] and references therein). Since in an 
environment of particles (holes) it can only move to the left (right) it always finds 
positions with positive density gradient (0. . . 021. . . 1). Its dynamics was defined in this 
way to localize the position of a shock (that is defined as a sudden change in the density 
approaching two different values to the left and right). Through the exact solution one 
was then able to calculate the density profile seen from the second-class particle. This 
has been considered as a limiting case of a shock profile (since on a ring the density is 
constant, the density seen from the second-class particle far to the left and right is the 
same). This form of a shock can also be found in one phase of the defect ASEP, namely 
for 1 — a > p > (3 (which does not occur for a = 1) and in the open-boundary ASEP 
along the second-order transition line pQ. Originally it described shocks with the same 
profile in the ASEP on the infinite line [T3] . 

In the continuous-time limit the critical density fl86l) becomes p c = (3 and for the 
occupations around the defect one obtains the well-known results from the defect ASEP. 
One finds in phase 1 (p > (3) that p_ = p and p + = 1 — (1 — (3)(1 — p) and for the 
defect velocity 1 — (3 — p. In phase 2 (p < (3) one has p_ = p 2 / [3, p + = 1 — (1 — p) 2 and 
v — 1 — 2p. The density profile has been calculated in [12] and one can take the results 
for the present process. In phase 1, where the defect behaves like a normal particle the 
density profile in front decreases exponentially to its bulk value p and behind the defect 
the density is constant. In phase 2 it behaves like a second-class particle and the density 
in front (behind) is increased (decreased) and reaches its asymptotic value algebraically. 

For parallel dynamics the ASEP with a single defect has not a natural equivalence, 
since the evolution of configurations in which the pattern 120 occur are not well-defined 
since under parallel dynamics 1 and 2 can not move to the right at the same time. 
However the process considered in this article solves this situation. 120 corresponds to 
10100. This moves into 01100 (210) at rate (3(1 - p) (12 exchange), into 10001 (102) 
at rate p{l — (3)) (20 exchange) and into 01001 (201) at rate pf3 (12 exchange, then 10 
exchange). We expect that the profile (in phase 2) of the current process plays a similar 
role for parallel dynamics. 
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7. Conclusions 

We shortly reviewed the solutions of the asymmetric simple exclusion process (ASEP) pQ 
and showed that the steady state weights for open boundaries can be written simply as a 
product of a scalar (pair-factorized) factor containing the nearest-neighbor correlations 
of the parallel update and a matrix-product state. In the second part we investigated 
a process on a ring in which particles have a maximum velocity 2, i.e. they can either 
move one site if they have exactly one empty site in front, or they can move two sites 
if they have more free sites in front This dynamics leads to an extinction of odd-valued 
gaps between consecutive particles. For overall odd number of holes in the system 
there remains with time exactly one excess hole that leads to a natural parallel defect 
dynamics. The presence of the excess hole in the odd case leads again to a product of 
a scalar pair-factorized and a matrix-product state which we assume to be generic for 
this type of driven- diffusive systems. The model exhibits a first-order phase transition 
separating two regimes with different defect velocities that we calculate exactly. As a 
step towards the calculation of the phase-dependent density profiles we obtain the exact 
expressions for the occupations behind and in front of the defect. 

We have shown how the model can be mapped onto the defect ASEP in the random- 
sequential limit. The ASEP with a second-class particle (being a special case) [13] turned 
out to have a density profile around the defect that could be considered as the limiting 
case (equal densities to the left and right) of a microscopic shock profile. It showed that 
this microscopic shape occurred in different ASEP contexts pQ. We expect that the 
present process plays a similar role for parallel dynamics. It seems to be the simplest 
process on a ring with one particle-species and conservative totally asymmetric dynamics 
with short-range interactions leading to a non-trivial steady state with phase transition. 

We finally want to point out that the process can be considered as a special case 
of a simple traffic model [16] in which particles can also move one site with a different 
probability a if it would be possible to move two sites (100 — > 010). A work on this in 
general ergodic model is in progress [TOl [15] . 
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Appendix A. Derivation of the normalization-generating function 

For the S n occurring in (j5J) we find 

S n = (l-p)(W 1 \G n {x,y)\V 1 )+py(W 1 \E 1 G n ^(x,y)\V 1 ), (A.l) 

where So = (1 — p)(l — 13). The functions G n obey the following recursions G n = 
C\G n -i + pxyKGn-i and respectively G n = G n _iCi + pxyG n -%K, with GLi := and 
G = 1, so that Gi = C 1; G 2 = C\ + pxyK, G 3 (x, y) = Cf + pxy[C x K + KCi) and so 
on. Here one has if = (1 — p){D\ + E\ + p) and C\ = xD\ + yE\. Special cases of the 
G{n) occurred in the open boundary case [5] with x — y — 1. It turns out that for x,y 
general it is difficult to work directly with G n (x, y). Instead we consider the generating 
function, which can be written as 

oo oo 

F(x, y, A) = A " G n(x, y) = £ A n + pxyXKf . (A.2) 

n=0 n=0 

The term under the sum is C% + pxyXK = (x + pxyX)D\ + (y + pxyX)E\ + p 2 xy\. It is 
very convenient to transform the matrices. Define primed matrices through 



One can check that these primed matrices indeed fulfill D' l E[ = (1 — p)(D[ + E[+p) = 
(1 — p)K'. In this notation C\ + pxyXK becomes finally 



C\ + pxyXK = y (x + pxy\)(y + pxyX)K' + lu, (A. 5) 

with 



to = uj(x, y, X) = \J (x + pxyX) (y + pxyX) (p — 2) + (1 — + y + 2pxyX) + p 2 xy X.(A.G) 
After execution of the sum in (1A.2|) the result can be written as 

1 / X\/(x + pxyX)(y + pxyX) \ 

F(x,y,X) = -^—[l-^ V \ 1 • ( A ' ? ) 

1 — cjA \ 1 — UM / 

One further needs an expression for the action of D[ and on the boundary vectors. 
For powers D'i\V) one gets for example 

A calculation adaptable from the defect-ASEP in [TTJ then yields for S(X) = 
Z^\ n S n = (W 1 \(l-p + py\E 1 )F\V 1 ): 



S= l^p + Pl 



I f3(l-p)y y+pxyXj 



1+pyX (3 



1-7 i — / y+pxyx \ 



x+pxy\ J 



(A.9) 
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with a function 7 to be determined from 



7(1-7) X^J(x + pxy\)(y + pxyX) 



(A.10) 



1 — p(l — 7) 1 — uX 



The singularity of flA.9j) closest to the origin is in phase 1 the pole at 7 



1 + (p — [3)/ (3/(1 — p) J (x + pxyX) /(y + pxy\) and in phase 2 a square root singularity 



in 7 resulting from (1A.10I) . 
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